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We propose the concept of a space mission to probe the so called flyby anomaly, an 
unexpected velocity change experienced by some deep-space probes using earth gravity 
assists. The key feature of this proposal is the use of GNSS systems to obtain an increased 
accuracy in the tracking of the approaching spacecraft, mainly near the perigee. Two 
low-cost options are also discussed to further test this anomaly: an add-on to an existing 
spacecraft and a dedicated mission. 

1. Introduction The flyby anomaly 

During the past couple of decades, a few deep-space probes that used an Earth flyby 
have apparently displayed an unexpected velocity change after their gravitational 
assist. This has become known as the flyby anomaly. 
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The effect in question was detected in the residuals of the analysis performed 
on the Doppler and ranging data, which showed the impossibility of fitting the tra- 
jectory with a single hyperbolic arc, but allowed for a separate fit of the inward 
and outward paths. The additional observed velocity shift is highly localized at the 
perigee, where tracking through the Deep Space Network (DSN) is not available 
(with an approximate four hours gap) . The spatial resolution of the available recon- 
structions, resulting form the 10 s interval tracking, does not allow for an accurate 
characterization of the effect, so that no corresponding acceleration profile exists. 
Only the variation of the probes' velocity (vis-a-vis kinetic energy) is known. 

This flyby anomaly has so far been observed in the Galileo, NEAR, Rosetta, 
and Cassini missions . A summary of the Earth flybys observed since the 1990s 
is shown in Table [if 21 A detailed discussion of the two Galileo (1990 and 1992) 
and the NEAR (1998) gravity assists is available in the literature^-! This includes 
an analysis of the three earliest flybys where the anomaly was observed, with an 
account of the accelerations generated by different known effects, in an attempt 
to single out possible error sources. An estimated average acceleration associated 
with the flyby anomaly of the order of 10~ 4 m/s 2 is measured against the Earth 
oblateness, other Solar System bodies, relativistic corrections, atmospheric drag, 
Earth albedo and infrared emissions, ocean tides, solar pressure, ei<Pl, 



Table 1 . Summary of orbital parameters from Earth flybys during the last couple 
of decades. 



Mission 


Date 


e 


Perigee 






AVoo/'Uoo 








(km) 


(km/s) 


(mm/s) 


(io- 6 ) 


Galileo 


1990 


2.47 


959.9 


8.949 


3.92 ±0.08 


0.438 


Galileo 


1992 


3.32 


303.1 


8.877 


-4.6 ± 1 


-0.518 


NEAR 


1998 


1.81 


538.8 


6.851 


13.46 ±0.13 


1.96 


Cassini 


1999 


5.8 


1173 


16.01 


-2±1 


-0.125 


Rosetta 


2005 


1.327 


1954 


3.863 


1.80 ±0.05 


0.466 


MESSENGER 


2005 




2347 


4.056 


0.02 ±0.01 


0.0049 


Rosetta 


2007 




~ 2400 




~ 




Rosetta 


2009 




2481 




~ 





Subsequently this discussion was extended to other possible error sources, com- 
paring this 10 -4 m/s 2 figure with several additional unaccounted acceleration 
sources. These include the atmosphere, ocean tides, solid tides, spacecraft charging, 
magnetic moments, Earth albedo, solar wind and spin-rotation coupling. It is con- 
cluded that all of the considered effects are several orders of magnitude below the 
flyby anomaly^. 

A quick overview of the magni tude s of all effects discussed in the two previous 
paragraphs is compiled in Table [2?^ one sees that all listed effects (except the 
Earth oblateness) are orders of magnitude smaller than the required value. This 
raises the issue of possible errors in the gravitational model of the Earth. However, 
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attempts to solve the flyby problem by changing the related second dynamic form 
factor J2 have yielded unreasonable solutions, and are unable to account for all 
aybysP. 

Table 2. List of orders of magnitude of possible error 
sources during Earth flybys. 



Effect Order of Magnitude 

(m/s 2 ) 



Earth oblateness 


10" 


-2 


Other Solar System bodies 


10" 


-5 


Rclativistic effects 


10- 


-7 


Atmospheric drag 


10- 


-7 


Ocean and Earth tides 


10" 


-7 


Solar pressure 


10- 


-7 


Earth infrared 


10- 


-7 


Spacecraft charge 


10- 


-8 


Earth albedo 


10" 


-9 


Solar wind 


10" 


-9 


Magnetic moment 


10- 


15 



An empirical formula to fit the flyby relative velocity change has been proposed 
by Anderson et a as a function of the declinations of the incoming and outgoing 
asymptotic velocity vectors, Si and S , respectively 

-jp = K (cos Si - cos 5 ), (1) 

where the constant K is expressed in terms of the Earth's rotation velocity uje, its 
radius Re and the speed of light c as 

1ui e R e 

K = — — . (2 

c 

This identification is suggestive, as it evokes the general form of the outer metric 
due to a rotating body-^, 

rfs 2 = f 1 + 2^-P\ (c dtf -fl- (dr 2 + r 2 dtf), (3) 



with 



$0 V 1 fu: P .R e s 1 



c 2 c 2 2 V c 



(4) 



where d£l 2 = d9 2 + sin 2 9d(f> 2 , Vq is the value of the Newtonian potential V(r) at 
the equator, uj e is Earth's rotational velocity and R e is its radius. Following this 
reasoning, and given the strong latitude dependence of Eq. (TTJ, this expression 



January 4, 2012 1:43 WSPC/INSTRUCTION FILE Galileo"FlybylJMP8 



4 O. Bertolami, F. Francisco, P.J.S. Gil & J. Paramos 

appears to indicate that the Earth's rotation may be generating a much larger 
effect than the frame dragging predicted by General Relativity. This, however, is in 
contradiction with the recent measurements of this effect performed by the Gravity 
Probe B probe^, which orbits the Earth at a height of ~ 600 km, well within the 
onset zone of the reported flyby anomaly. 

2. Effect on GNSS systems 

In order to discuss the possible use of current and future GNSS constellations to 
probe this putative flyby anomaly, one should first evaluate to what extent it can 
affect their individual elements. These follow approximately circular Medium Earth 
Orbits (MEO), at a height of about ~ 20000 km; since the anomalous velocity 
change is only observed before and after flybys occurring at much smaller heights 
(of the order of 1000 km), one may empirically dismiss any effect. 

One could sharpen the above argument, even though a full analysis is impossible 
due to the lack of spatial resolution and consequent inability to fully characterize 
the spatial dependence of the reported anomaly^. Notwithstanding, one takes as 
relevant figure of merit the anomalous acceleration a ~ 10~ 4 m/s 2 , which may be 
assumed constant in the absence of further data. In this case it has been shown that 
no constant acceleration greater than 10~ 9 m/s 2 can affect the GNSS constellation, 
since it would have otherwise been detected^. 

Thus, one concludes that the flyby anomaly, if real, must be due to a strongly 
decaying force, which should drop by four orders of magnitude with a modest (about 
fourfold) increase in distance, from r = Re + h ~ 7000 km to r ~ 27000 km. 
As a result, one may safely assume that the GNSS constellation is fundamentally 
unaffected by this putative anomaly, and may be thus employed to track probes 
performing gravity assists at the relevant region h ~ 1000 km. 

3. GNSS spacecraft tracking 

The tracking of spacecraft through GNSS systems is already commercially available 
(e.g. EADS-Astrium's MosaicP, NASA PiVoT^). These systems are typically used 
to follow satellites in low earth orbit (LEO), at altitudes below those of the GNSS 
satellites (h < /ignss ~ 20000 km), where the GNSS signal is stronger. Nevertheless, 
the Equator-S mission can receive front lobe signal from GPS satellites at an altitude 
of 61000 kirM Furthermore, it is worth exploring the possibility of using the side 
and back lobes of the GPS 

signaL J10111] 

to establish non-line of sight tracking and 
avoid the shading of the Earth. Clearly, the build up of more constellations and the 
use of multi-GNSS receivers, able to work simultaneously with different systems, 
will increase the accuracy of above-MEO satellite tracking in the coming years. 

The accuracy of GNSS spacecraft tracking is, understandably, better for lower 
orbits; however, it should be noted that during the apogee of highly elliptical or- 
bits (HEO), the velocity is, of course, much slower than close to perigee. This 
allows for the construction of a good orbital solution, despite the decreased signal 
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coverage^ES As a result, the position and velocit y accuracies for different types 
of orbit are somewhat similar, as depicted in Table 

Recall that there is no full characterization of the anomalies during the flyby, and 
these are detected from the mismatch between the expected and observed velocities 
after gravitational assist. As stated before, this is due to the inability of the DSN 
to track the spacecraft trajectories very close to the atmosphere, during a ~ 4 h 
gap. Regarding the possibility of using the GNSS in this region, Fig. Q] shows that, 
although the velocity error is maximum close to perigee, this peak is very localized: 
from a baseline of ~ 20 mm/s during the remaining orbit, it peaks briefly at ~ 
100 mm/s (during the first perigee approach), and converges towards ~ 50 mm/s 
in the subsequent perigee passings. By plotting the aforementioned gap, one sees 
that accuracies of ~ 20 mm/s are attainable during approximately half of this time 
interval. 

For the study of the flyby anomaly, one would be interested in a high velocity 
accuracy, at least of the same order of magnitude as the observed Avoo ~ 1 mm/s. 
The currently available systems provide around 20 mm/s, which is clearly insuffi- 
cient for such a study. However, the presented accuracies are related to real-time 
orbit solutions, which is unnecessary for the purpose of this study, and can undoubt- 
edly be improved if offline processing is used, alongside other weak signal tracking 

[To] 

strategies^ 21 . This, together with the increasing numbers of elements of the available 
(and upcoming) GNSS, lead us to conclude that it is indeed feasible to use the latter 
to test the flyby anomaly, if not with the current capability, then in the near future. 



4. Options for probing the flyby anomaly 

We consider two options to test the flyby anomaly: an add-on to an existing mission 
on a Highly Elliptic Orbit (HEO), or a dedicated low-cost mission in either HEO 
or a hyperbolic trajectory. 
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Fig. 1. Velocity error of multi-GNSS tracking of HEO spacecraft. Boxes (centered on perigee 
with 4 hour width) signal the gap in DNS coverage; the horizontal line corresponds to a 20 mm/s 
accuracy (adapted from Ref. (12)). 
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Table 3. Typical accuracies expected from GNSS satellite tracking 
systems for LEO, MEO, Geosynchronous Earth Orbit (GEO) and 
HEO. 



Orbit 


Apogee 
(km) 


Position Accuracy 
(m) 


Velocity Accuracy 
(mm/s) 


LEdni 


200 to 2000 


10 


10 


MECQI 


2000 to GEO 


30 


20 


GEdH 


35786 


150 


20 


HECUSI 


> 35786 


100 


20 



In the first option, the choice would be to piggyback a multi-GNSS receiver in 
an existing space mission. Sinc e these receivers are relatively compact and with 
reduced power consumption^^, the host mission could be a small low-cost one. At 
perigee, a highly elliptical trajectory would present a comparable (although smaller) 
velocity and height as the reported anomalous gravitational assists, with the added 
benefit of allowing for repeated experiments. 

One can take as an example of a suitable mission the cancelled Inner Magne- 
tosphere Explorer (IMEX) mission of the NASA University Explorer programme, 
with a mass of only 160 kg and a total budget of 15 Mp3 in 2 000. The IMEX 
probe was to be launched as a secondary payload on a Titan IV launcher, but was 
cancelled due to cost overrun. It would have followed a HEO, as summarized in 
Table |31 which would provide a "flyby" velocity at perigee of about 10 km/s, close 
to the reported anomalous flybys. 

The more ambitious option of a dedicated mission naturally has a number of 
advantages over the former, the main of which is the choice of orbit that can closely 
mimic a gravity assist, including an hyperbolic one. However, as discussed above, 
a closed orbit of sufficiently high ellipticity would provide for multiple flybys, in- 
creasing the quality of the obtained data and allowing for a better characterization 
of the anomaly. The HEO would also allow to ascertain if the flyby anomaly is 
exclusively linked to hyperbolic orbits. Also, possible error sources such as aerody- 
namic and thermal effects close to perigee could be more closely controlled with 
a dedicated mission. For instance, the spacecraft could be enclosed in a spherical 
radio-transparent body, so to simplify modelling and reduce directional effects. If a 
spin is given, any accidental anisotropics would be averaged out, yielding a much 
cleaner testbed for the desired experiment. 

This mission would require a micro-satellite with a mass under 100 kg and a bud- 
get cap similar to the IMEX mission. This upper bound is rather straightforward 
to argue by comparison. Firstly, no additional spending is anticipated, due to the 
simplified spherical design over the more complex IMEX probe. Secondly, the sci- 
entific instrumentation found in the latter would be replaced by just a multi-GNSS 
receiver, thus lowering the total cost. More ambitiously, an added accelerometer 
could provide for a cost-effective independent measure of the acceleration profile, 
with a modest addition to the mass budget (such as the ~ 3 kg /iSTAR instrument 
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considered in the Outer Solar System and Odissey mission proposal d^^ . 

Following in the purpose of this paper, to present the feasibility of using GNSS 
to probe the fiyby anomaly, this estimate illustrates the low cost of a dedicated 
mission for that purpose. Nevertheless, the actual cost could, in principle, be some- 
what smaller than around 15 M$, the IMEX cost estimate, not only due to the 
inherently simpler design and instrumentation, but also due to the ongoing trend of 
decreased micro-satellite costs, reflecting advances in miniaturization, lower power 
consumption and improved industrial processed^. 

Table 4. Orbital parameters of IMEX's Highly Elliptical Orbit 
and similar hyperbolic flyby. 





IMEX 


Similar Hyperbolic Orbit 


Perigee altitude 


349 km 


349 km 


Apogee altitude 


35800 km 




Velocity at perigee 


10.1 km/s 


11 km/s 


Eccentricity 


0.7248 


1.04 


Orbital period 


10.5 h 





5. Conclusions 

In this work the use of the Galileo system or the GNSS to study the flyby anomaly 
is proposed. One finds that most of the available studies dealing with the tracking 
of spacecraft in real time have an insufficient velocity and position accuracy to 
detect this discrepancy. However, since this real time accuracy is only one order 
of magnitude above the required one (in particular, ~ 10 mm/s vs. ~ 1 mm/s in 
velocity), it is reasonable to expect that this situation could improve in the short- 
term. A thorough exploitation of available resources could lead to a suitable tracking 
of spacecraft with greater accuracy, by abandoning real time solutions, and resorting 
instead to offline processing, use of side and back lobe tracking, amongst other weak 
signal tracking strategies. Crucially, the use of several GNSS at once should lead to 
an increased coverage of the different geometries. 

Thus, it can be safely stated that there is no fundamental issue preventing the use 
of GNSS tracking to study the reported flyby anomaly. Naturally, this availability 
is not sufficient, as only spacecraft equipped with a (multi-)GNSS receiver would 
allow for such a study. In this work, we have shown that a mission of this kind 
could be easily deployed, either as an add-on package to an existing platform with 
the required highly elliptical orbit, or through a dedicated mission. While the first 
scenario would provide a cheap solution, we argue that a dedicated mission could 
be envisaged with a higher scientific payoff, while maintaining an overall low-cost 
approach. 

Regardless of the actual origin of the flyby anomaly (unaccounted conventional 
effect, precision glitch or, more appealingly, new physics), we believe that our pro- 
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posal offers a low-cost opportunity for displaying some of the scientific possibilities 
opened by the GNSS era. 
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